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Abstract 

We present a theoretical study of the interaction of vacancies with a tilt grain boundary in 
aluminum based on the density functional theory. The grain boundary volume expansion and 
vacancy induced contraction are calculated and compared for the nearest-neighbor atoms, and the 
former is found to be smaller than the latter. The formation energy of a vacancy placed at various 
layers in the grain boundary has been calculated and we find that the grain boundary does not 
always act as sinks for vacancies. In fact, it costs more energy to form a vacancy at the boundary 
plane than in bulk, although the rest of the grain boundary region does attract vacancies. The 
microscopic mechanisms of grain boundary sliding and migration are investigated thoroughly with 
and without a vacancy. We find that although the vacancy can hinder the grain boundary motion 
by tripling the energy barrier of sliding and migration, it can not inhibit or even delay the migration 
process. The vacancy placed at the first layer from the interface is found to be trapped at the layer 
and not able to follow the migrating interface. 



I. INTRODUCTION 



Vacancies are the simplest lattice defects and play an important role for various prop- 
erties of materials, such as kinetic, thermodynamic, electrical, optical and even mechanical 
properties, to name a few. Similarly, grain boundaries as planar defects, are also responsible 
for the various material properties. [1], |2|] While there have been extensive theoretical studies 
aimed at either vacancies or grain boundaries individually, much less effort has been devoted 
to the understanding of the collective behavior of grain boundaries and vacancies, especially 
from parameter-free ab initio studies. For example, ab initio calculations have been used 
with great success to study the vacancy formation energy and the vacancy-induced relax- 
ation in different metals, p, |J the vacancy diffusion energy barrier and the annihilation 
process in semiconductors, || || [?]] etc. However, all these calculations are confined to either 
bulk systems or surfaces. On the other hand, grain boundaries (pure or with impurities) 
have attracted a great deal of attention recently. For example, ab initio calculations have 
been performed for the grain boundary energies and the atomic structures of the pure Ell 
tilt and E3 twin grain boundaries in aluminum ||, and the effect of impurity on the grain 
boundary cohesion in Ni 3 Al. || However, these studies do not involve vacancies in the grain 
boundaries and examine only the properties of the grain boundaries themselves. Therefore, 
one of the goals in this paper is to investigate the combined effect of grain boundaries and 
vacancies in aluminum using ab initio calculations and to explore the possible interactions 
between the point and planar defects. 

While the equilibrium structures of grain boundaries have been extensively studied, only 
very limited research has focused on the atomistic simulation of grain boundary sliding and 



migration. [ 10|, |TT|, |12|, |13| The former process involves the relative displacement of the two 
constituent grains in a direction parallel to the boundary interface, while the latter involves 
the motion of the interface in a direction perpendicular to the boundary plane. Grain 
boundary sliding/migration is considered to be one of the principal mechanisms of plastic 
deformation of polycrystalline materials at intermediate-to-high temperatures (above 0.6T m , 
where T m is the melting point) and it contributes to creep deformation and intergranular 
fracture, which can in turn be significantly affected by the presence of vacancies. [0, [14| These 
two modes of grain boundary motion can be correlated from a simple geometrical analysis. 
T5| Very recently, Molteni et al. carried out ab initio studies of the atomic structure and the 
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energy barriers for grain boundary sliding in germanium and aluminum. [ 16 , [I7jj Employing 



molecular dynamics simulations in conjunction with the Embedded Atom Method, Chandra 
and Dang have shown that the coupling between sliding and migration depends strongly on 
the external conditions that initiate the grain boundary sliding. [T2[ 



In this paper we have studied the interactions of vacancies with the tilt £5 (210) [001] 
grain boundary in aluminum using state-of-the-art ah initio total-energy calculations based 
on the pseudopotential plane wave method. First, we have calculated the vacancy forma- 
tion energies and atomic structures for a vacancy placed at various positions in the grain 
boundary region. These results reveal that the grain boundary may act as sink or source of 
vacancies depending on the location of the vacancy. Second, in order to understand the ef- 
fect of vacancies on the grain boundary mobility, we have investigated the effect of vacancies 
on the grain boundary sliding and migration behavior. Aluminum is selected in this study 
because it is a prototypical simple metal of importance in industrial applications and some 
of the conclusions we draw from the present study may be extended to other fee metals as 
well. 

The remainder of this paper is organized as follows: Sec. II describes briefly the compu- 
tational techniques and the supercell used for the grain boundary calculations. In Sec. Ill 
we present the results for the atomic structures and energetics of the pure grain boundary. 
Sec. IV contains the results with vacancies introduced at various atomic sites of the grain 
boundary. The vacancy formation energy at different sites is calculated and compared to 
that of the bulk system. The atomic relaxation due to the grain boundary expansion and 
the vacancy-induced contraction will be examined and compared. In Sec. V we present 
the energy profile and the atomic structures for the grain boundary sliding and migration 
process with and without a vacancy. We find that the vacancy can significantly increase 
the grain boundary sliding and migration energy barrier, and therefore hinder the grain 
boundary motion. In order to shed light into the microscopic physics of the combined grain 
boundary sliding/ migration motion, we have examined in detail some critical atomic config- 
urations during the sliding/migration process. In Sec. VI a brief summary and statement 
of conclusions are presented. 

II. MODEL AND COMPUTATION 

The electronic structure of the grain boundary is calculated by means of the pseudopo- 
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tential plane wave method based on the density functional theory and the local density 
approximation. We employed the exchange and correlation potential of Ceperly and Alder 



as parametrized by Perdew and Zunger [ 18 1 and the norm-conserving nonlocal pseudopoten- 
tial of Bachelet, Hamann and Schliiter. [T9|] The Kohn-Sham wave functions are expanded 
in plane waves with a kinetic energy cutoff of 14 Ry. The atomic structures are consid- 
ered fully relaxed when the Hellmann-Feynman forces on each atom are smaller than 0.001 
Ry/au. For the Brillouin-zone sampling, a fc-point grid consisting of (2,7,1) divisions along 



the reciprocal-lattice directions is used according to the Monkhorst-Pack scheme. P0[| Con- 
vergence tests have been performed both for the number of A;-points as well as the number 
of plane waves. The calculated equilibrium lattice constant (ao), bulk modulus and cohe- 
sive energy are 3.93 A, 80.59 GPa and 3.26 eV/atom, respectively, in good agreement with 
the corresponding experimental room-temperature values of 4.05 A, 76.93 GPa and 3.39 
eV/atom. 

In this work we consider the symmetric £5 (210) [001] tilt grain boundary in the 
coincident-site-lattice (CSL) notation. In the CSL notation, the grain boundary is formed 
through rotation of the two fee grains about the [001] tilt axis by 36.9°. The boundary is 
oriented along the (210) plane and the inverse density of CSL sites, £ = 5, characterizes 
the unit cell volume of the CSL superlattice. The E5 [001] (210) grain boundary of Al is 
simulated by the orthorhombic supercell shown in Fig. 1. The supercell contains 60 atoms 
distributed in 30 (210) layers parallel to the interface. The large and small circles represent 
the atoms on the (001) and (002) planes, respectively. The dimensions of the supercell along 
the [120], [001] and [210] directions are 8.78 Ax 3.93 Ax 26.35 A, respectively. Because of 
the periodic boundary conditions, the supercell contains two equivalent grain boundaries. 
However, as will be shown in Sec. Ill, the separation between the two grain boundaries is 
large enough to eliminate the interaction between them. For the modeling of vacancies in the 
grain boundary, a single vacancy is placed at all possible independent lattice sites labeled by 
an integer n = 0-7, which represents the nth nearest-neighbor layer from the grain boundary 
plane (the number of independent lattice sites is reduced due to the symmetry between the 
(001) and (002) planes and that between the two grain boundaries). In addition to the grain 
boundary calculation, we have also computed the vacancy formation energy for the Al bulk 
by employing a simple cubic supercell with 108 atomic sites to model an isolated vacancy. 
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III. PURE GRAIN BOUNDARY 



A pure grain boundary in the following is referred to as the grain boundary without 
vacancies. The atomic structure of the pure grain boundary before and after relaxation is 
shown in Fig. 2. The open and filled circles represent the atomic positions before and after 
the relaxation and the large and small circles correspond to the atoms on the (001) and (002) 
planes, respectively. The grain boundary volume expansion is clearly seen near the interface 
and one local measure for such an expansion is the relative normal displacement of the two 
atomic planes closest to the interface. The present calculation yields a local expansion of 
0.09 do, where do is the bulk Al lattice constant. This local expansion of the Al atoms in the 
pure Al £5 grain boundary is smaller compared to the corresponding expansion of the Al 
atoms in the M3AI alloy for the same grain boundary. || This is consistent with the fact that 
there is charge transfer from the Al atoms to the nearest-neighbor Ni atoms in M3AI, which 
in turn reduces the screening charge between the nearest Al-Al ions across the interface, 
leading to a stronger electrostatic repulsion between the Al ions in M3AI grain boundary. In 
Fig. 3 we present the relative variation (strain) of the interlayer spacing as a function of the 
layer away from the grain boundary plane. The strain has a symmetric oscillatory profile 
that reaches its maximum at the first layer and decays into the bulk. The nature of this 
strain profile can be traced to the Friedel oscillations of the electron distribution which drive 



the ions to relax in such an oscillatory pattern. [^TJ The negligible strain for the layers from 
6 to 10 clearly indicates that the atomic relaxation is localized only within 5 layers from the 
interface. Therefore, the supercell with 15 layers between the two interfaces is large enough 
to capture the overall properties of the grain boundary. 

It is interesting to notice that even after the relaxation, the distance between the nearest- 
neighbor atoms across the interface is 2.44 A, about 20% smaller than the corresponding 
distance of 2.78 A in the bulk system. This result invalidates the common misconception 
that the grain boundary is associated with open space. On the other hand, as will be shown 
in Sec. IV, vacancies can be formed relatively easier in most of the grain boundary region 
compared to the bulk, giving rise to more empty space in the region. The grain boundary 
energy can be determined from the difference of the energy of a unit cell containing the grain 
boundary and a unit cell containing an equal number of atoms in the bulk environment, 
divided by the grain boundary area. We find a grain boundary energy of 502 mJ/m 2 for the 
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E5 grain boundary, which is larger than the corresponding EAM value of 351 mJ/m 2 .f22| 
This result is consistent with the finding that, overall, the EAM potential tends to give 
smaller grain boundary energies compared to ab initio calculations. || 



IV. GRAIN BOUNDARY WITH A VACANCY 

The relaxed atomic structures of the grain boundary before and after the introduction of 
a vacancy are shown in Fig 4. (a), (b), (c), corresponding to the vacancy placed at the zero-th 
(interface), first and seventh layers from the interface. The large and small circles correspond 
to atoms on the (001) and (002) planes, respectively. The open (filled) circles represent the 
relaxed atomic positions before (after) the introduction of the vacancy, and the arrows show 
the direction of the atomic displacements. In order to gain insight into the interplay between 
the grain boundary and the vacancies, we have calculated the displacements due to the grain 
boundary expansion (u^s) and that due to the vacancy contraction (uy). Both vlgb and 
uy are computed for the nearest-neighbor atoms (first shell) from the defects. For the case 
of the vacancy at the interface (Fig. 4(a)), although the atomic relaxation is small and 
localized, the vacancy-induced contraction (uy = -0.12 A) is much larger than the grain 
boundary expansion (ugb = +0.01 A). When the vacancy is placed at the first layer, the 
atomic relaxation is more dramatic (Fig. 4(b)). There are seven atomic layers involved in the 
relaxation triggered by the vacancy and the first shell displacements are uy = -0.53 A while 
ugb= +0.34 A, respectively. In Fig. 4(c) we show the atomic structure with the vacancy 
placed at the seventh layer (bulk like). We find uy = -0.06 A, which is about 2% of the 
nearest-neighbor distance of the bulk lattice and is along the [110] direction. These results 
are in excellent agreement with the calculations by Turner et al. |J for a vacancy in Al bulk 
both for the magnitude and the direction of uy. In other words, the atomic structure for 
a vacancy at the seventh layer from the interface completely reproduces the characteristics 
of a vacancy in the bulk. Along with the fact that the grain boundary expansion ugb at 
the seventh layer is zero as discussed earlier, we conclude that the selected supercell for the 
grain boundary calculations is large enough to yield reliable results. It is also interesting 
to note that as far as the nearest-neighbor shell is concerned, the vacancy-induced atomic 
relaxation is more pronounced than that due to the grain boundary. 

In order to clarify whether grain boundaries act as sinks or sources for vacancies, we 
have calculated the vacancy formation energy for various layers from the grain boundary 
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interface. The results are presented in Fig. 5, where the straight line denotes the vacancy- 
formation energy in bulk. The striking feature of the energy curve is the large reduction of 
the vacancy formation energy at the first layer, which strongly suggests the preference of 
forming a vacancy at the first layer compared to bulk. Although the figure also shows that 
some other more distant layers in the grain boundary are also energetically favorable for the 
vacancy formation, one can clearly see that the interface itself is not favorable for forming 
a vacancy. In fact, there is more than 0.1 eV energy penalty for creating a vacancy at the 
interface compared to bulk. Note that the vacancy formation energy at the fifth layer from 
the interface is close to the bulk value, showing the convergence of the energy with respect 
to the supercell size. In conclusion, the ab initio calculations support the traditional view 
that grain boundaries act as pathways for vacancy diffusion and that vacancies prefer to 
segregate to grain boundaries over bulk at low stress. However, we should also point out 
that one out of five lattice sites in the grain boundary region turns out to be energetically 
unfavorable for vacancies, namely those sites on the grain boundary plane. 

V. GRAIN BOUNDARY SLIDING AND MIGRATION: EFFECT OF VACANCIES 

In this section we present results for the grain boundary sliding and migration, and the 
effect of vacancies on the grain boundary motion. 

In our ab initio modeling, the grain boundary sliding is simulated quasi-statically, namely 
the top grain is rigidly shifted over the bottom by a series of small specified distances along 
the interface. The sliding distance is described in percentage of a cs i, where a cs i is the lattice 
parameter of the CSL cell along [120]. Clearly sliding by a cs i will bring the grain boundary 
back to its original configuration due to the imposed periodic boundary conditions. For 
each rigid shift, atomic relaxation is performed to locate the closest local energy minimum 
configuration. In order to understand the effect of vacancies on the grain boundary sliding 
and migration process, we have calculated the grain boundary energy profile as a function of 
sliding distance for the pure grain boundary and the grain boundary with a vacancy placed 
at the first layer from the interface, corresponding to the lowest vacancy formation energy 
at the grain boundary. 

The grain boundary energy profile during the sliding process is shown in Fig. 6, where 
the filled circles represent the pure grain boundary and the open circles represent the grain 
boundary with the vacancy, respectively. The labels (a)-(f) on both curves indicate the 
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special atomic configurations which will be explored in detail below. For the pure grain 
boundary, there are two energy peaks at 30% and 80% a cs i, the first much broader than 
the latter; and one energy valley at 55% a cs i. It is important to note that the vacancy 
triples the energy barrier for sliding and migration, which demonstrates the great hindering 
effect of the vacancy on the interface mobility. However, the detailed analysis of the atomic 
structures below shows that the vacancy can not prevent or even delay the occurrence of 
the grain boundary migration. 

The collective motion of atoms during the grain boundary sliding and migration is essen- 
tial to understand this seemingly trivial, but in fact complicated process. To shed light on 
the related physical process, we next examine the atomic configurations associated with the 
special points (a)-(f) in the energy curves of the pure grain boundary in Fig. 6. The relaxed 
atomic structures corresponding to (a)-(f) are shown clockwise in Fig. 7, with the filled and 
open circles representing the atoms on the (001) and (002) planes, respectively. Note that 
more atoms are displayed in these figures than the actual atoms used in the supercell, which 
is depicted by the solid frame. In order to monitor when does grain boundary migration take 
place, one needs to locate the position of the interface during the sliding process. However, 
since the grain boundary is not symmetric any more during the sliding process, it is not al- 
ways obvious to find the new position of the migrating interface. To overcome this difficulty, 
we propose a simple scheme which allows us to identify the moving interface by following the 
evolution of the structural unit shown in the figures. The characteristics of this unit is that 
it consists of two sub-units that cross the interface with a common vertex on the interface. 
These sub-units are the two-dimensional unit cells of the distorted fee lattice. By locating 
such a unit and the common vertex during the sliding process, one can identify the position 
of the grain boundary plane. Applying this scheme, we find that the migration of the pure 
grain boundary starts at 35% a cs i and ends at 55% a cs i. After the completion of the first 
migration, the grain boundary has to return to its original position at 100% a cs i to complete 
the second migration. For the grain boundary with the vacancy, the migration starts at the 
exact same sliding distance, 35% a cs i, but ends at about 62% a cs i. The second migration 
due to the periodic boundary conditions is also observed as in the pure grain boundary. 

Fig. 7(a) shows the initial atomic structure for the pure grain boundary before sliding. 
The interface is indicated by a solid line. The angle a (ZABC) of the structural unit is slightly 
distorted from its value of 90° in the perfect fee lattice. By examining the various atomic 
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configurations before the first migration (35% a cs /), we find that they are all similar, being 
different mostly along the sliding direction. More specifically, the angle a of all configurations 
turns out to be obtuse and one example of such a configuration is shown in Fig. 7(b). When 
the interface starts migrating, however, the grain boundary experiences an abrupt overall 
structural transformation, as illustrated in Fig. 7(c). The interface has clearly moved from 
its original position (dotted line) to the current solid line position, and the angle a has 
changed from obtuse to acute. This one-layer migration reflects the geometrical necessity 



for the coupling between migration and sliding. |T|, [15| After the onset of the migration, 
the atomic structures remains akin to each other during the sliding until the migration is 
fully completed, in which, as shown in Fig. 7(d), the grain boundary recovers its original 
symmetric structure. Note, that the angle a now changes from an acute angle to a right 
angle. The dotted and solid lines are the same as in Fig. 7(c). Because of the periodic 
boundary conditions, the grain boundary has to return to its original position during the 
remaining sliding process. The critical configuration it has to cross is shown in Fig. 7(e), 
where the grain boundary experiences the highest energy. The intriguing feature of this 
configuration is that no boundary plane can be identified, while there is a "mirror" plane 
(dashed line) across which the atoms on the (001) plane are symmetrically situated to those 
on the (002) plane. This half-way plane reflects the transient position of the migrating 
interface and therefore the energy associated with this configuration gives the energy barrier 
for the grain boundary migration. Finally, Fig. 7(f) represents the completion of the sliding 
process where the grain boundary returns to its initial symmetric configuration. The dotted 
line indicates the position of the interface after the first migration and the solid line shows 
the current position (also the initial position before the sliding) of the interface. 

For the case of the vacancy placed at the first layer from the grain boundary, the sliding 
energy profile becomes more irregular, with multiple peaks and valleys (Fig. 6). The abun- 
dance of the fine grain boundary structures reflects the complexity of the system, so that we 
shall concentrate only on the most important events. Fig. 8(a) shows the atomic structure 
before sliding and the vacancy is denoted by the letter V. The solid line indicates the current 
position of the interface. It is found that this initial state no longer corresponds to the lowest 
energy, which is reached after a sliding displacement of 2% a cs j. There is an energy cross-over 
at 10% & cs i, above which the grain boundary energy with the vacancy becomes larger than 
the corresponding value of the pure grain boundary, until the next cross-over occurs at about 
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75% a cs /. Although the energy value does not return to zero as in the case of the pure grain 
boundary when the migration is completed, careful examination of the atomic structures in 
Fig. 8 indicates the occurrence of the grain boundary migration. Just like the pure grain 
boundary, the atomic structures remain similar to each other during the sliding process until 
the first migration takes place at about 35%. One such configuration is shown in Fig. 8(b), 
where the angle a is obtuse. The highest energy configuration is reached when the interface 
migrates one atomic layer upward, shown in Fig. 8(c), in which the dotted and solid lines 
represent the original and current positions of the interface, respectively. As in the case of 
the pure grain boundary, the atomic structures at the onset of the migration differ signif- 
icantly from those before the migration and the atoms near the interface become severely 
distorted from their positions in the underlying fee lattice. The grain boundary completes 
its first migration at about 62% a cs /, where the symmetric structure is recovered (Fig. 8(d)). 
It is important to notice, however, that this symmetric structure with the vacancy (labeled 
by V) at the interface is different from the initial state where the vacancy was at the first 
layer above the interface. Therefore, the energy difference between the initial state (a) and 
the current state (d) reflects simply the difference in the vacancy formation energy between 
the vacancy placed at the interface and the first layer, respectively. This indicates that the 
vacancy is trapped in the first layer and not able to follow the migrating interface. The grain 
boundary develops another peculiar configuration at 92% & cs i, just before the completion of 
the second migration. A symmetric pattern forms relative to the previous interface (dashed 
line) and the free volume associated with the vacancy disappears. The off-stoichiometry due 
to the vacancy, however, is recovered at the symmetry plane (dashed line). The solid line 
represents the final (also initial) position of the interface during the entire sliding process. 
The atomic structure of the grain boundary at 100% & cs i is shown in Fig. 8(f) with the letter 
V indicating the vacancy. The dotted (solid) lines represent the grain boundary position 
after the first (second) migration. 

VI. CONCLUSION 

We have studied the interaction between vacancies and the S5 tilt grain boundary in Al, 
using ab initio pseudopotential calculations, which yield reliable results for both the atomic 
structure and energetics compared to empirical atomistic approaches. The nearest-neighbor 
distance across the interface is actually smaller than that in the bulk even after atomic 
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relaxation. We have calculated the grain boundary volume expansion and the vacancy- 
induced contraction and found that the latter is larger than the former for the nearest- 
neighbor shell. We have also calculated the vacancy formation energy for a vacancy placed 
at different layers from the interface. We find that the grain boundary does not always act 
as sinks for vacancies. While most of the grain boundary sites have lower vacancy formation 
energy compared to bulk, those at the interface have higher formation energy. We have 
investigated the combined grain boundary sliding and migration process with and without a 
vacancy. We have introduced a simple scheme to identify the migrating interface during the 
sliding process. We have found that although the vacancy can hinder the grain boundary 
motion by tripling the energy barrier for sliding and migration, it can not inhibit or delay 
the occurrence of migration. In fact, the grain boundary starts the migration at 35% a cs / 
regardless if there is a vacancy or not. We have analyzed the atomic configurations during 
the sliding process in order to gain better insight on this complicated process. Interestingly, 
the vacancy is found to be trapped at the first layer and not able to follow the migrating 
interface. 
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FIG. 1: Supercell for the £5 (210) [001] tilt grain boundary in Al viewed along the [001] direction. 
The large and small circles represent atoms on the (001) and (002) planes, respectively. The 
integer (n = 0-7) represents the nth nearest-neighbor layer from the grain boundary interface that 
a vacancy is placed. The grain boundary plane is indicated by the solid line. 

FIG. 2: Unrelaxed (open circles) and relaxed(filled circles) atomic structure for the pure £5 grain 
boundary. Large and small circles represent atoms on the (001) and (002) planes, respectively. 

FIG. 4: Relaxed atomic structures for the grain boundary before (open circles) and after (filled 
circles) the introduction of a vacancy at the interface (a), the first layer (b) and the seventh layer 
(c) from the interface, respectively. The large and small circles represent atoms on the (001) and 
(002) planes, respectively. The arrow shows the direction of the atomic displacement due to the 
vacancy. 

FIG. 5: Vacancy formation energy as a function of the number of the layer away from the interface 
on which a vacancy is placed. The vacancy formation energy in bulk Al is denoted by the solid 
line. 

FIG. 6: Relative grain boundary energy as a function of sliding distance expressed as percentage 
of the lattice parameter, a cs i, of the CSL cell along [120]. The filled and open circles correspond to 
the energies of the pure grain boundary and the grain boundary with a vacancy at the first layer, 
respectively. The labels from (a) to (f) represent the special atomic configurations displayed in 
Figs. 7 and 8. 

FIG. 7: Evolution of the atomic structure of the pure grain boundary during the sliding process. 
The snap shots of the grain boundary structures correspond to the special points (a)-(f) on the 
energy curve for the pure grain boundary in Fig. 6. Filled and open circles represent atoms on 
the (001) and (002) planes, respectively, and the arrow indicates the direction of grain boundary 
migration. 

FIG. 3: Interlayer strain (percentage) as a function of the number of layers away from the interface. 
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FIG. 8: Evolution of the atomic structure of the grain boundary during the sliding process with 
a vacancy placed at the first layer above the interface. The snap shots of the grain boundary 
structures correspond to the special points (a)-(f) on the energy curve for the grain boundary 
with the vacancy in Fig. 6. Filled and open circles represent atoms on (001) and (002) planes, 
respectively and the arrow indicates the direction of grain boundary migration. 
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